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We present a simple theory for structural transitions in 
displacive ferroelectrics of the perovskite type. In our the- 
ory, the competition between the elastic energy cost for the 
displacement of the homopolar ion from the centrosymmetric 
position, and the energy gain due to a ferroelectric order- 
ing of the dipoles formed by the ionic displacements, leads 
naturally to a first-order transition from a paraelectric to a 
ferroelectric phase. This transition takes place at a certain 
temperature T C (L) as the temperature is decreased and, at a 
certain size L C (T) as the size of the system is increased. The 
transition temperature T C (L) is suppressed as the sample size 
is reduced, and vanishes for samples below a certain size. For 
order-disorder ferroelectrics, our theory shows that the sup- 
pression of T c by a reduction in system size is not appreciable, 
a result that is borne out by experiments. 

PACS Numbers: 68.35 Rh, 77.84 -s, 77.84 -Dy 



I. INTRODUCTION 

Though the experimental study of finite pSize effects 
in ferroelectric materials has a long historjoB, the rapid 
development of advanced synthetic techniques has now 
made it possible to study different compounds in the form 
of phase-pure, ultrafine particles with a narrow size dis- 
tribution. There is also a strong motivation for studying 
size-limited ferroelectric systems in view of their current 
and potential applications as sensors, memory elements, 
nano-robotic and micro-electromechanical devicesa. 

A ferroelectric is termed 'displacive' when the elemen- 
tary dipoles strictly vanish in the paraelectric phase, and 
'order-disorder' when they are non-vanishing but ther- 
mally average out to zero in the paraelectric phase. In 
recent years, the displacive ferroelectric transitions in 
PbTiC>3 and BaTiC>3 nanoparticles have been studied in 
detail. The nature of the size effects observed in both 
systems is essentially similar. A decrease in the parti- 
cle size causes monotonic reductions in (a) the transition 
temperature T c , and (b) the tetragonal distortion of the 
unit cell which characterizes the ferroelectric phase. So, 
at a low enough particle size, the lattice tends to assume 
the high temperature, high symmetry, cubic paraelectric 
structure. p 

Uchino et alB have found that T c (in degrees Centi- 
grade) falls with particle size L (nm) following the rela- 
tion 

T e (£) = T e (oo) - f , (1.1) 

where T c (oo) = 128 (500), B = 700 (588.5), and L c = 
110 (12.6) for BaTiOa (PbTi0 3 ). 



In a recent study of nanocrystalline PbTiOsJjsing di- 
electric, thermal and structural measurements^], it was 
established that with decreasing particle size: (1) there 
is a monotonic decrease in the T c , (2) the value of the 
peak dielectric constant (e max ) decreases, (3) the ferro- 
electric transition becomes increasingly diffuse, and (4) 
the crystallographic unit cell tends towards higher sym- 
metry (c/a — » 1). Qualitatively similar results have-been 
obtained for PbZrOa , a displacive antiferroelectriaH. 

It is also instructive to compare finite size effects in 
displacive and order-disorder ferroelectrics. In electri- 
cally insulated samples of sub- micron KH2PO4 , a typical 
order-disorder system, the depolarization field appeared 
to prevent th£_stabilization of ferroelectric ordering below 
a critical sizaiil. However, a later study of NaN02 showed 
clearly that there was no suppression of T c down to 5 nm 
in samples suspended in either electrically insulated or 
conducting mediaO. It is clear that the size driven effect 
is much weaker for order-disorder ferroelectrics. This is 
a feature that our theory also captures (Section [v]). 

In recent years, there have been a few attempts 
to understand theoretically the nature of size effects 
in ferroelectrics. Using the pheupmenological Landau- 
Devonshire theory, Zhong et al&3 have shown that the 
ferroelectric T c should decrease with decreasing size, ul- 
timately leading to a size-driven phase transition fropj 
the ferroelectric to the paraelectric phase. Shih et alM 
have considered the effect of incorporating the depolar- 
ization energy in the Landau free energy density. 

In this paper, we restrict ourselves mainly to the de- 
scription of quasi-free ferroelectric nanoparticles. The 
system is assumed to consist of loosely aggregated, un- 
damped particles which are not electrically isolated. Un- 
der such circumstances, we can neglect the effects of ex- 
ternal strain and depolarization. This approximates the 
experimental situation considered by Chattopadhyay et 
al.U. Specifically, we have selected PbTi03 as the model 
system, but the results should apply to other displacive- 
type systems as well. 

PbTiC>3 is a classical displacive ferroelectric with a 
tetragonal perovskite structure (a = 0.3899nm, c = 
0.4153nm). At room temperature, the Ti and O ions are 
displaced with respect to the Pb ions, parallel to the polar 
axis with displacements d-pi = -018nm, cIq = .047nmt3. 
On increasing the temperature, the tetragonal ferroelec- 
tric undergoes a first order transition to a cubic paraelec- 
tric phase (a — c — 0.396nm) at T c = 768K. 

The paper is organised as follows. In section |l| we 
present our model and discuss the qualitative physics of 
the phase transition. Section [II gives the details of the 
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calculation, and our results are presented in section IV. a \ _ y^fi^, £ 2 ■ ifc 4 £ 4 + acr 2 — ba f 2 ) 



In the picture adopted in the present paper, the tran- 
sition takes place as a result of competition between an 
ordering tendency of dipoles in adjacent cells and the 
elastic energy cost associated with the displacement of 
the atoms. Our main results are the following. We offer 
a simple physical explanation (see Fig.|J) for the first- 
order transition in displacive systems driven by tempera- 
ture and system size. We calculate the strain of the unit 
cell and make comparison with experiments. While the 
temperature dependence of the strain (Figj^) is not in 
very good agreement with experiments, the size depen- 
dence we obtain is in excellent agreement (Fig.|5|). Wc 
show that for order-disorder ferroelectrics, size effects are 
highly suppressed (section.^). We present the phase di- 
agram for both kinds of ferroelectrics in the temperature 
- size plane in Fig|| The paper ends with a brief dis- 
cussion of certain features not captured quantitatively in 
our theory, and prospects for future work (section |VT]) . 



II. THE MODEL AND THE PHYSICAL PICTURE 

A dipole moment is created in each unit cell by the 
motion of the ionized atom at the centre of the unit cell 
to an off-centre position. This atom is known as the ho- 
mopolar atom. In PbTiC>3 and BaTiC>3 for instance, the 
Ti atom plays this role. The displacement is accompa- 
nied by a distortion of the cubic unit cell to a tetragonal 
one with sides a x & x c. The order parameter measured 
in the experimentsa, a = cj a — 1, obtained from powder 
x-ray scattering, is a measure of the strain or distortion 
of the unit cell. 

Clearly there is an energy cost at each site associated 
with the displacements £j of the homopolar atom. In 
spite of this cost, the system may find it profitable to un- 
dergo a distortion accompanied by an off-center motion of 
the homopolar ion if there is a negative energy contribu- 
tion from the interaction energy of dipoles when they are 
aligned. We describe this energy phenomenologically by 
an interaction —J£,i£,j between nearest neighbours. We 
refer to this term as the 'Ising' term (see below). For 
ferroelectrics, J is positive, whereas a negative J will de- 
scribe antiferroelectrics. In the present work we focus 
mainly on ferroelectrics (see however section VI). 

The effective Hamiltonian for the problem is thus the 
sum of an elastic part and an Ising part. The form of 
the elastic part may be deduced from simple symmetry 
considerations. The cost of a displacement ^ of the ho- 
mopolar atom at the i th unit cell and an associated strain 
<7j(= c/a — 1) of the unit cell can be written as a power- 
series expansion in £j and <7i, at each site. Only even 
powers of £j are allowed in the expansion since the cost 
cannot depend on which way the atom moves. There is 
no such restriction for the strain cr, which is a scalar. We 
thus arrive at the Hamiltonian 



(2.1) 



There is no term proportional to <Zj alone, since when 
£i = for all i, the system should be in equilibrium with 

(I, = 0. 

An additional consideration has gone into the trun- 
cation of the power series. It is expected on symmetry 
grounds that the thermally averaged value (cr), of the 
strain will be an even function of the thermally averaged 
displacement (£), and for small distortions, (cr) ~ (£ 2 ). 
This is borne out extremely well from experimental data 
on PbTiOa (Fig. |l]) . From the data we find the empirical 
relation 
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where A = 3.24 for PbTiOa . The quadratic relation 
holds good for other oxides as wellliU, with different values 
of A. 

We may conclude from the above discussion that £ 2 
and a are of the same order of smallness. The power 
series expansion in (2.1) thus retains all terms upto 0(a 2 ) 
or C>(£ 4 ). 



Starting with the Hamiltonian (2.1), we can easily in- 
tegrate out the strain field cr, to arrive at an effective free 
energy in terms of the displacements, 



i i 
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(2.3) 



where A2 = &2 an d A4 
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FIG. 1. Experimental values of the strain (a = - — 1) plot- 
ted as a function of d/a, the displacement |£| of the homopolar 
atom in units of lattice spacingli3. The points correspond to 
measurements at differ ent temperatures below T c , and the 
line is the best fit (Eq. 
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FIG. 2. Schematic plots of the different contributions to 
the free energy ( ^.3|) . The elastic part (dashed line) has a 
minimum at the origin. The second minimum is produced by 
the Ising part (dotted line) which is constant up to a certain 
displacement and begins to decrease after that. This results 
in the full free energy with two minima (solid line). 



For very large L, z — > 6 and the system, let us say, is fer- 
roelectric provided T < T c . If we reduce L, this reduces 
z(L), and the magnitude of the Ising part decreases as 
a result. This makes the ferroelectric phase unstable at 
small size and the system becomes a paraelectric. 

In particular, at zero temperature, when there are no 
thermal fluctuations, cti is the same at all sites, and the 
free energy is simply 

f(d) = ±(\ 2 -Jz(L))d 2 + ±\ i d\ (2.7) 

which describes a second order transition driven by size 
when 

Jz{L c ) = A 2 (2.8) 
at a certain critical size L c . 

III. CALCULATIONS 



To get a physical picture of the temperature driven 
transition in bulk samples, we introduce variables on and 
di corresponding to the direction and magnitude of the 
displacement, 

6 = dion, (2.4) 

where on (= ±1) is an Ising variable. With this change of 
vari able s, the ferroelectric part of the Hamilfconian (See 
Eq.( |3.l|) ) is like a compressible Ising terrnifj, since the 
coupling ~ Jd 2 depends on the displacement. At a tem- 
perature T, the Ising part cannot be ordered (and hence 
the system cannot be a ferroelectric) if the coupling is less 
than T. This means that the system can be a ferroelectric 
only ifd 2 is greater than ~ T / J . Whether it actually is a 
ferroelectric or a paraelectric depends on energetics. This 
explains why we expect a first-order transition driven by 
temperature. 

To capture this quantitatively, we may perform a par- 
tial trace over the a degree of freedom and evaluate a 
free energy functional of the distortions di alone, given 
by 

e- F ^=Tr ai e- H ^ d l (2.5) 

The global minimum of F(d), determines the thermody- 
namic phase. We shall see in the following section that 
this procedure automatically leads to a first order tran- 
sition. 

The size induced transition is also easily understood 
within this picture. For a system with N — (L + l) 3 
sites, the elastic energy is ~ N and the Ising energy is 
~ Nz(L) where z(L) is the average coordination number 
of the cubic lattice of linear size L, given by 

4L) = {L l 1)3 [6(i - l) 3 + 30(L - l) 2 

+ 48(L - 1) + 24] . (2.6) 



In this section we present a mean-field theory for our 
model. The effective Hamiltonian ( |2.3| ) expressed in 
terms of di and a.i is 

F(d,a)= 1 -X 2 J2d 2 + 1 -X i J2 d t 

i i 

— J didjOLiCtj. (3-1) 

<ij> 

Recall that di is the magnitude of the displacement of the 

th 

homopolar ion in the i unit cell and on is its sign. We 
now perform a trace over the a variables at the mean- 
field level, to get 

i i 

+ — J mirrijdidj 

<ij> 

. , JdiJ^ 4 m^di 
-T^>2cosh( ^ 3 J ). (3.2) 

i 

The thermal average m, = (q,;) is determined by the self 
consistency equation 

Jdi V . m-jdj 
m, = tanh( ^ J ). (3.3) 

In this paper, we resort to the approximation di = d, 
i.e., we assume that the homopolar atom displaces by 
the same amount in all the unit cells. Therefore, the free 
energy per unit cell f(di) — F(di)/N is given by (see Eq. 

CD) 

f(d) = e(d)+I(d), (3.4) 
where we have separated the elastic part 
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e(d) = ^A 2 d 2 + ^d 4 , 



(3.5) 



and the Ising part 

n 1 V 1 ^ Jd 2 V 

i(d) = jd 2 ~Y. m ^ - ^E Tln2cosh ( — t 



<ij> 



(3.6) 



For a bulk system (L — > oo), the Eq.(3.3) gives a uniform 
solution mi = m. We thus have 

7(d) = - Jzm 2 d 2 - Tin 2 cosh(^^). (3.7) 

where z(— 6) is the coordination number of the cubic 
lattice. 

The schematic plots of e(d),I(d) and f(d) are shown 
in Fig. |^. It can be seen that the elastic part is an in- 
creasing function of the displacement with a minimum 
at d = 0, whereas the Ising part is constant, equal to 
—Tin 2, for d < d c = ^JT/Jz and decreases (i.e. in- 
creases in magnitude, being negative always) for larger 
d. Indeed, for d < d c , to = (see Eq. (3.2)) and 
therefore 1(d) = —Tin 2. As a result of this behaviour 
of 1(d), the free energy develops a second minimum at 
some d = do, which is found by solving the equation 
df/dd — 0. For large d, 1(d) ~ — d 2 , so that the 
quartic term in e(d) dominates for large d, resulting in a 
stable free energy function as shown in Fig. |2| It is this 
two-minima structure of the free energy that gives rise 
to a first-order transition between the ferroelectric (the 
d = do minimum) and paraelectric (the d = minimum) 
phases. 



Minimising the free-energy in Eq.(3.4) with respect to 
d, we get 



\ 2 d + A 4 d 3 - Jzm d = 0, 

corresponding to the two minima, one at d 
other at 



d 



' Jzm 2 — A 2 



A 4 



(3.8) 
and the 



(3.9) 



At the first-order transition, the two minima coexist, i.e., 
/(0) = /(do), or 



1 



1 



1 



-A 2 do + 7A 4 do + - Jzm 2 d^ - T c lncosh( 



Jzm 2 d\ 



= 0. 



(3.10) 



We need three conditions to fix the three parame ters 
of o ur th eory A 2 ,A4, and J. The last two equations ( 3.9 
and 3.10) supply two of these conditions; in these we use 
the experimental values of T c and d . The third condition 
is the equation A 2 = Jz(L c ) derived at the end of the last 
section. 



From experiments on PbTiC>3 , we have do =j-0.299A, 
T c — 768K, and from the work of Zhong et alL3, L c ~ 
10. Using these in the three equations above, we obtain 
J = 4.345909 x 10 4 KA~ 2 , A 2 = 2.3718 x 10 5 Kl~ 2 , and 
A 4 = 2.6529 x 10 5 Kl~ 4 . 

Now it is not surprising that since Jzd 2 /T c ~ 30.27, 
we get m ~ 1 for all the sites (see Eq.(3.3)). This in fact 
turns out t o be the ca se als o for a finite lattice, so that the 
equations (|3.9D and (3.1C) can be used even in this case 



by replacing z by z(L). The results from this completely 
agree with the results of solving Eq.(3.3) explicitly for 
a finite lattice with open boundary conditions, which we 
have performed numerically. 

To obtain the phase diagram, w e cal culate do(T) for 
each L from equations (3.£) and ( 3.10 ) (with z(L) re- 
placing z). 



IV. RESULTS 

Fig.|| shows our result for the strain as a function of 
temperature for the bulk case (solid line); the experimen- 
tal datatZl are also shown (crosses). The lack of quantita- 
tive agreement is ascribed to the fact that in our mean- 
field theory, the ferroelectric is completely saturated, i.e., 
to ~ 1 at all T < T c . 
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FIG. 3. The strain plotted as a function of temperature 
(T) for bulk samples. The crosses are the experimental points 
and the solid line is from the theory. The dotted line is ob- 
tained by assuming a weak linear temperature dependence for 
A2, and the dashed line is by assuming a temperature depen- 
dent A 4 (Fig.§). 

We can make up for this disagreement in two ways. 
One possibility is to assume a weak temperature depen- 
dence for A 2 , of the form AT + B\ For B = 5.45 J, A = 
5.5 x 1CT 4 J, A 4 = 1.23JA~ 2 and J = 1.65 x 10 5 Kl- 2 
(T measured in degrees Kelvin) this leads to the dotted 
line shown in Fig.[| We find that while this makes the 
order parameter temperature dependent, it also makes 
it J-independent and therefore size- independent, though 
the size-driven transition still remains. 

Alternatively, we can assume a temperature depen- 
dence for A4 which can be calculated by requiring that 
the order parameter agree with experiment. This turns 
out to be a rather strong temperature dependence for A4 
which we show in Fig.|4[ Our subsequent results are ob- 
tained by assuming that the parameters A2 and A4 are 
independent of temperature. 

We have shown in Fig|| a plot of strain calculated as 
a function of system size at a temperature q£ 300 K, 
along with the corresponding experimental datau. In this 
graph, we have normalized the calculated strain in a way 
so as to match the experimental value for the bulk sys- 
tem at this temperature. We obtain excellent agreement 
with experiment for the size dependence of the strain. 




FIG. 4. The ratio R = A 4 (T) /A 4 (T C ) plotted as a function 
of temperature. This dependence leads to the dashed line 
shown in Fig|| 
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FIG. 5. The strain in the tetragonal phase (- — 1) as a 
function of reduced system size at room temperature. The 
solid line shows the strain calculated from our theory (see 
text). The crosses are experimental pointsB. 

Fig.p shows our phase diagram in the temperature - 
size plane, obtained by calculating strain. The phase 
boundary for displacive systems (solid line) is a line of 
first-order transitions for all T > 0, and the size-driven 
transition at T = is continuous, as noted above. It 
can be seen from the phase diagram that our theory pre- 
dicts a suppression of T c as the system size is reduced, in 
agreement with experiments on displacive systems and 
also other theoretical studies on these systems. 



FIG. 6. The phase diagram of ferroelectric oxides in the 
temperature-size plane obtained within our theory. We have 
plotted the reduced temperature t = T C (L) /T c (oo) against the 
reduced particle size. The solid line is for displacive systems 
and is a first-order line separating the tetragonal ferroelectric 
phase from the cubic paraelectric phase. It can be seen that 
there is an appreciable suppression of T c as the system size 
is reduced, and T c is zero below a certain size. The dotted 
line is for order-disorder systems and is a line of second-order 
transitions separating the low-temperature ferroelectric phase 
from the high temperature paraelectric phase. The transition 
temperature T c remains practically equal to its bulk value 
down to very small system sizes, and the size effects are highly 
suppressed in this case. In particular, the system remains in 
the ordered phase below T c /2 for all sizes. 



V. ORDER-DISORDER FERROELECTRICS 



For order-disorder ferroelectrics, the local potential 
experienced by the homopolar atom has a double well 
structures, i.e., there is a local maximum at d = and 
two minima at d = ±do- This corresponds to A2 < 
within our approach. 

In this case, it is easy to see that size effects are highly 
suppressed. Firstly, we note that the Ising term can never 
lead to a minimum at d = at any finite temperature. 
Therefore an increase in temperature can destroy the fer- 
roelectric order via a continuous transition. T he d istor- 
tion continues to have the value given in Eq.(3.9), but 
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rn = (a) vanishes when Jz(L)d^/T = 1. The last con- 
dition can be met either by changing T or by changing 
system size. For the bulk system (z = 6), this is satisfied 
for T = T c . We can therefore write z{L) /6 ~ T/T c , where 
L is the size at which there is a size-driven transition at 
temperature T. Note that when expressed in terms of 
the reduced transition temperature t = T C (L)/T c (oo) the 
above equation (which describes the boundary between 
the ferroelectric and paraelectric phases) is independent 
of our model parameters. If T/T c is about 0.9 (say), 
the n we have to go to such small systems as L ~ 10 (see 
Eq.( [2.6| )) to see size effects. 

Typically, the lattice spacings in these systems are ~ 
4 A, so we have to get samples of size ~ 4nm, to have 
observable size effects, even at so high a temperature as 
0.9T C . At lower temperatures the critical sizes are even 
smaller. In particular, for T < 0.5T C , we do not expect 
any size-driven transition, since_s(L) cannot be less than 
3. The experiment on NaNOjlj that we referred to in 
Sec.[| reports that T c is not suppressed down to sizes of 5 
nm. 

We have shown the phase diagram for order-disorder 
ferroelectrics also in Fig.^, and the difference between 
the two kinds of ferroelectrics is very instructive. We at- 
tribute the suppression of size effects in the case of order- 
disorder ferroelectrics to the absence of a structural tran- 
sition accompanying the ferroelectric transition. Indeed, 
at any temperature, one has to go to much smaller sam- 
ples to meet the requirement J z{L)d^/T = 1 (which is 
the transition condition for order-disorder systems) than 
the co-existence condition Eq.( [3.10 ) (which is the transi- 
tion condition for displacive systems). 



VI. DISCUSSIONS 

Although our model gives a useful physical picture of 
both temperature and size driven transitions, it is rather 
too simple to provide perfect quantitative agreement. We 
have made several simplifying assumptions, the most sig- 
nificant being that the displacement has the same mag- 
nitude throughout the sample. This is not necessarily 
true, especially for small systems where the displacement 
at the surface could be different from that it is in the 
bulk. This, and the subsequent neglect of thermal fluc- 
tuations, turn out to be rather drastic approximations. 
The strongly first-order character of the transition (i.e. 
the very weak, almost nonexistent temperature depen- 
dence of the order parameter) is probably an artifact of 
our mean-field approximation. 

An interesting question, which cannot be answered 
within the framework of our theory, is whether the tran- 
sition remains first-order beyond mean-field theory. Fur- 
ther, since our theory is an effective single-site theory, it 
cannot capture capture the phenomenon of mode soften- 
ing that is observed in these systemscj. 

We do not have the effect of depolarization field in our 



theory. The depolarization field, under certain circum- 
stances, caftJead to better quantitative agreement with 
experimented. However, in the experimental situation 
that we are describingp, the effect of the depolarization 
field should be negligible. While we have not made ex- 
plicit the role of long range interactions between dipoles, 
this will simply lead to a reparametrisation of J within 
the simple version of our theory presented in this paper. 

The mean-field theory that we have attempted on the 
Ising part of the Hamiltonian consists in reducing the lat- 
tice problem to an effective single-site problem. This led 
to a very weakly temperature order parameter as we have 
seen. We have also performed, calculations based on the 
Bethe-Peierls approximation^, which takes into account 
short range correlations between dipoles. However, this 
did not lead to any -better temperature dependence for 
the order parameter^. 

There are certain qualitative merits of our theory that 
are worth emphasizing. For example, our theory can be 
easily extended to include antiferroelectric oxides, sim- 
ply by having a negative J. The ordered phase will have 
a nonzero value for an antiferroelectric order parameter 
such as sublattice magnetization as- The molecular field 
at any site i, due to the nearest neighbours, will be oppo- 
site in sign to the order parameter at that site, so that the 
product of the molecular field and J will have the same 
sign as for the ferroelectric case. The theory as worked 
out above in terms of displacements will go through with- 
out any further change, and the results will be identical 
to that for ferroelectrics, with J replaced by \ J\. Recent 
experimental studies of size-driven transitions in antifer- 
roelectric materials!] report results very similar to those 
as ferroelectrics, confirming our e xpec tation. 

We can explain the form of Eq.(l.l) for the size-driven 
transitions quite simply in o ur th eory. For the size-driven 
transition, we get from Eq.( 3.10 ), upon replacing A2, A4, 
and J by their values determined from bulk data, an 
equation of the form lncosh(z(L)/T' c ) = z(L)/T c — A, for 
some constant A. It is cle ar that this equation will have 
a soluti on o f the form fll.lD for large L since z(L) is given 
by Eq.( |2.6|) . The actual numbers appearing in Eq.([L~l|) 
will of course depend on the system we address. 

Our theory has three fitting pa rame ters, whereas the 
phenomenological Landau theorieslHML.3 have about twice 
as many. This circumstance is clearly because our the- 
ory is based on a microscopic model which identifies the 
different contributions to the system free energy as being 
due to the elastic and ferroelectric parts. As a result, we 
are also able to describe the structural transition in the 
displacive systems in addition to the ferroelectric transi- 
tion, to which it is related. The chief qualitative merit 
of our theory is that the first-order transition appears in 
a very natural way, with minimal assumptions about the 
form of the elastic part of the Hamiltonian. With the 
lowest order nonlinear elastic term (i.e. fourth order), 
we are able to describe the first-order transition. This is 
an essential difference with earlier approaches based on 
a Landau expansion of the free energy. 
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Finally, we have shown that our theory can describe 
size-driven transitions in order-disorder systems as well 
as displacive ones. Our simple estimate for the critical 
size at which size effects become important in these sys- 
tems is in excellent agreement with data on NaNC^lii 
Our approach is superior to the Landau theory for size- 
driven transitional^ since size effects in order-disorder 
systems can be addressed in our framework. 

In conclusion, we have presented a simple unified pic- 
ture for structural transitions in displacive and order- 
disorder ferroelectric oxides, which captures the physics 
of both temperature-induced transitions in the bulk and 
size induced transitions in nanoparticles. The size effects 
are understood in a simple way as being the result of 
smaller coordination number near the surface of the par- 
ticle. Our mean-field approximation gives qualitatively 
good results, except for the temperature dependence of 
the order parameter. Future work will center around in- 
corporating the effect of thermal fluctuations and getting 
closer agreement with experiments. 
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